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can be taken over, some of the scaling relations between the ageing exponents are 
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fluctuation-dissipation ratio. The form of the scaling function of the two-time response 
function is compared with the prediction of the theory of local scale-invariance. 
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1. Introduction 

Symmetry principles have been extremely useful for the understanding of complex many- 
body systems, where the interactions between the degrees of freedom are sufficiently 
strong as to render perturbative methods inapplicables. Here we are interested in 
the slow non-equilibrium dynamics shown by many-body systems which are rapidly 
brought out from some initial state ('quenched') to a region in phase-space where 
either the equilibrium state naturally generates a slow dynamics (this is for example 
realized for systems at a critical point) or else into a coexistence region dominated by 
several equivalent stationary states. One of the essential features of such systems is that 
their properties depend on their 'age', that is the time elapsed since the quench. Of 
course, any biological system ages, but there is also 'physical ageing' which arises even 
if the underlying microscopic dynamics is completely reversible. One might formally 
define ageing by this breaking of time-translation invariance, associated with a slow 
dynamics which generically leads to some form of dynamical scaling. Physical ageing 
was originally seen to occur in glassy systems [120] and has been used since prehistoric 
times by engineers in the processing of materials. Quite recently, it has been realized 
that very similar phenomena can also be found in simple magnets, without disorder nor 
frustrations. The study of these supposedly simpler systems may lead to conceptual 
insights which in turn could become also fruitful in more complex systems. The topic 
has been under intensive study, see [ISl E31 E21 [32l ESI EH EH EH [70] for reviews. 
The reversibility of the microdynamics in systems undergoing physical ageing means 
that their stationary states are equilibrium states. In numerical simulations of such 
systems this is realized by choosing the dynamics such as to satisfy detailed balance. 
In many systems undergoing physical ageing, detailed balance and consequently the 
relaxation towards equilibrium steady-states is taken for granted, as there are many 
textbook proofs [1231 1129j of detailed balance for closed, isolated systems. 

On the other hand, it has become increasingly clear from studies in anomalous 
chemical kinetics that several of the constitutive properties of ageing are naturally met 
in many situations. First, it is well-known that fluctuation effects may lead to slow, non- 
exponential relaxation in irreversible chemical reactions - not accounted for by mean-field 
schemes, see e.g. |117[ IM] for reviews and references therein. Second, it was understood 
more recently through the work of Oshanin and collaborators [971 l98l l6ll3T | 1 124(1125] that 
even for reversible reactions a slow, non-exponential relaxation may generically occur 
without the fine-tuning of parameters and furthermore, that the steady-states to which 
relaxation occurs depend on the kinetic coefficients and hence cannot be equilibrium 
states. Consequently, detailed balance cannot be valid in these systems. While in 
these studies long-range interactions (as they may naturally arise in reactions of large 
molecules or in studies of radiation damage [6]) play an essential role, detailed balance 
need not be satisfied even in kinetic systems with contact interactions. For example, 
in the system defined by the simultaneous reversible reactions 2A < — > 0, 2A < — > A, 
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A < — > with diffusive motion of single particles, detailed balance only holds is certain 
conditions on the reaction rates are met [I] . Since detailed balance is already found to 
be broken in very simple reactions such as 2 A < — > B [1241 1125] or A + B — > P,[3T], it 
is conceivable that the phenomenon might be much more common. Furthermore, since 
sometime a slow relaxation in kinetic models is brought into relationship with glassy 
dynamics [90] , it is of interest to investigate to what extent the three essential properties 
of physical ageing - slow dynamics, dynamical scaling and breaking of time-translation 
invariance - are actually realized in chemical kinetics. 

Therefore, we shall review here recent progress about ageing in many-body systems 
with a more general dynamics where detailed balance is no longer required to hold 
and therefore non-equilibrium steady-states may arise. Since there are as yet only few 
studies available on these systems, comparison with ageing in simple magnets should 
be a useful guide. For a similar reason, we shall investigate the behaviour of models 
whose physical origin is very different which should lead to some insight about generic 
properties. Because of the possibility of non-equilibrium steady-states, the systems 
under consideration here are closer to biological/chemical ageing than those considered 
up to now in studies of physical ageing. Remarkably, there is evidence that some 
dynamical symmetries recently discovered in physical ageing may also extend to this 
more general class of systems. 

We shall define the systems we want to study in the next section. Before we 
come to that, however, we shall briefly recall for reference some of the main results 
about the ageing of simple magnets. We assume throughout that the order-parameter 
is non- conserved by the dynamics and that the initial state is totally disordered, unless 
explicitly stated otherwise. Besides the breaking of time-translation invariance, ageing 
systems are often characterized by dynamical scaling. It has become common to 
study ageing behaviour through the two-time autocorrelation and (linear) autoresponse 
functions 



where <p(t, r) is the order-parameter at time t and location r and h(s, r) is the conjugate 
magnetic field at time s and location r. The scaling behaviour is expected to apply in 
the so-called ageing regime where t, s 3> t m icro 

and t S S> tmicroj 

where t m i cro is a 

microscopic time-scale. We illustrate this in figure Q] which shows the autocorrelation 



quenched from a fully disordered initial state to a temperature T < T c . When plotting 
the data over against t — s, we see that the data depend on both t — s and s, hence 
time-translation invariance is broken and the system ages. Further, with increasing 
values of the waiting time s, the system becomes 'stiffer' and a plateau close to the 

f Recall that the zero-temperature Glauber model can be mapped, via a duality transformation [118] 
or a similarity transformation to the kinetic model 2 A — ► with single-particle diffusion. 




(1) 
(2) 



function C(t,s) of the 3D Ising model with (non-conserve' 



d) Glauber dynamicgjl [19] 
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Figure 1. (a) Ageing and (b) dynamical scaling of the two-time autocorrelation 
function C(t, s) in the 3D Glauber-Ising model quenched to T = 3 < T c , for several 
values of the waiting time s [105] . 



equilibrium value C eq = develops when t — s is not too large before the correlations 
fall off rapidly when t — s —>■ oo. When replotting the same data over against t/s, a 
data collapse is found if s is large enough which is evidence for dynamical scaling. 

The distance of such systems from a global equilibrium state can be measured 
through the fluctuation-dissipation ratio, defined as [33] 

X(t,s):=TR(t,s)(^^y 1 (3) 

At equilibrium, X(t, s) = 1 from the fluctuation-dissipation theorem. One often 
considers the limit fluctuation-dissipation ratio := lim^oo (lirn t _ KX j X(t, s))jj For 
quenches to below T c , one usually has Xoo = but for critical quenches onto T = T c , 
it has been argued by Godreche and Luck that X^ should be a universal number [5T] . 
since it can be written as a ratio of two scaling amplitudes. This universality has been 
thoroughly confirmed for systems relaxing towards equilibrium steady-states, see [321124] 
for recent reviews. 

Furthermore, in writing eqs. ([1121 it was tacitly assumed that the scaling derives 
from the algebraic time-dependence of a single characteristic length-scale L{t) ~ t l / z 
which measures the linear size of correlated or ordered clusters and where z is the 
dynamic exponent. For a 2D Glauber-Ising model quenched to T = T c the growth of 
correlated clusters in illustrated in figure [2] where the black/white site represent the two 
states of the Ising spins. Then the above forms define the non-equilibrium exponents 
a and b and the scaling functions fc(y) and fniy). For large arguments y — > oo, one 
generically expects 

| The order of the limits is crucial, since Unit—too (lim^oo X(t, s)) = 1. 
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Figure 2. Snapshots from the 2D Glauber-Ising model quenched to T — T c from a 
disordered initial state at (left panel) t = 25 and (right panel) t — 275 MC time steps 
after the quench [105] . 



fc(y) ~ y- Xc/z , f R {y) ~ <r Afl/ * (4) 

where Ac and A#, respectively, are known as autocorrelation [161 [79] and autoresponse 
exponents [103J. While in non-disordered magnets with short-ranged initial conditions 
one usually has Ac = \r, this is not necessarily so if either of these conditions is 
relaxed. From a field-theoretical point of view it is known that for a non-conserved order- 
parameter the calculation of Xc,r requires an independent renormalization and hence 
one cannot expect to find a scaling relation between these and equilibrium exponents 
(including z) [81]. On the other hand, the values of the exponents a and b are known. For 
quenches to T = T c , the relevant length-scale is set by the time-dependent correlation 
length L(t) ~ ~ t x l z and this leads to a = b = (d — 2 + rj)/z, where r\ is a standard 
equilibrium exponent. For quenches into the ordered phase T < T c , one usually observes 
simple scaling of C(t,s) = fc(t/s), hence b = 0§| The value of a depends on whether 
the equilibrium correlator is short- or long-ranged, respectively. These may be referred 
to as classes S and L, respectively and one has, see e.g. [301 ESJ E3] 

{ e -kl/€ ( c i ass s \ 1/z 

\r\-( d -^) jdassL ° = { (d - 2 + rj)/z (5) 

Examples for short-ranged models (class S) include the Ising or Potts models in d > 1 
dimensions (and T < T c ), while all systems quenched to criticality, or the spherical 
model or the 2D XY model below the Kosterlitz-Thouless transition are examples for 
long-ranged systems (class L). 

In equilibrium critical phenomena, it is well-known that the standard scale- 
invariance can, under quite weak conditions, be extended to a conformal invariance. 
Roughly, a conformal transformation is a scale-transformation r i— > br with a space- 
dependent rescaling factor b = b(r) (such that angles are kept unchanged). In particular, 

§ This needs no longer be the case when the ageing close to a free surface is considered [10] . 
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in two dimensions conformal invariance allows to derive from the representation theory 
of the conformal (Virasoro) algebra the possible values of the critical exponents, to set up 
a list of possible universality classes, calculate explicitly all n-point correlation functions 
and so on [T6J [28] . One might wonder whether a similar extension might be possible at 
least in some instances of dynamical scaling and further ask whether response functions 
or correlation functions might be found from their covariance under some generalized 
dynamical scaling with a space-time- dependent rescaling factor b = b(t, r) [56], [58] ? We 
shall discuss the question here in the specific context of ageing and shall focus on what 
can be said about the scaling functions fc,R,(y) in a model- independent way. 

A useful starting point is to consider the symmetries of the free diffusion (or free 
Schrodinger) equation 

2Md t <l) = A0 (6) 

where A = V ■ V is the spatial laplacian and the 'mass' M. can be seen as a 
kinetic coefficient. Indeed, it was already shown by Lie more than a century ago 
that this equation has more symmetries than the trivial translation- and rotation- 
invariances. Consider the so-called Schrodinger- group defined through the space-time 
transformations 

. at + f3 . Rr + vt + a 

t^t = —-T I r^r' = — — , a8-p 1 = \ 7 

7^ + jt + d 

where a, /3, 7, 5, v, a are real (vector) parameters and R is a rotation matrix in d spatial 
dimensions. The group acts projectively on a solution of the diffusion equation through 
(t,r) g(t,r), <p^T g <f) 

(T 9 0)(t,r) = / 9 (^ 1 (t,r))0(^ 1 (t,r)) (8) 

where g is an element of the Schrodinger group and the companion function reads 
[951 [102] 



f g (t,r) = ( 1 t + 6)- d / 2 



cxp 



M jr 2 + 2 Rr ■ (7a — 8v) + 7a 2 — t5v 2 + 27a ■ v 
~2 -ft + 5 



(9) 



It is then natural to include also arbitrary phase-shifts of the wave function within 
the Schrodinger group Sch(d). In what follows, we denote by sd) d the Lie algebra of 
Sch(d). The Schrodinger group so defined is the largest group which maps any solution 
of the free Schrodinger equation (with M. fixed) onto another solution. This is easily 
seen in d = 1 by introducing the Schrodinger operator 

S := 2M y-! - Y 2 1/2 (10) 



The Schrodinger Lie algebra sct) 1 = (X_i i0 ,i, i , M ) is spanned by the infinitesimal 
generators of temporal and spatial translations (X_i,YL_i/2), Galilei-transformations 
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(X1/2), phase shifts (Mo), space-time dilatations with z = 2 (X ) and so-called special 
transformations (Xi). Explicitly, the generators read [56] 

X n = - t n+1 d t - ^Vrd r - n ^ + l) Mt n ~ l r 2 - -(n + l)t n 
2 4 2 

Y m = - t m+1/2 d r - (m + A^ m " 1/2 r (11) 

Here a; is the scaling dimension and A4 is the mass of the scaling operator <fi on which 
these generators act. The non- vanishing commutation relations are 

[-^ru^n'] = ( n ~ n ')X n+ri i , [X n , Y m ] = ^— — rrij Y n+m 

[X n ,M n >] = -n'M n+n , , [Y m ,Y m ,] = (m - m')M m+m , (12) 

The invariance of the diffusion equation under the action of scf^ is now seen from the 
following commutators which follow from the explicit form ffTTT) 

[S,X^\ = [S,Y ±1/2 ] = [S,M ] = 

[S,X ] = -S , [S,Xi] = -2tS - (2ar - 1)M (13) 

Therefore, for any solution <p of the Schrddinger equation S(j) = with scaling dimension 
x = 1/2, the infinitesimally transformed solution Xcp with X e sci) l also satisfies the 
Schrddinger equation SX(j) = [83], [95], [55] . For applications to ageing, we must consider 
to so-called ageing algebra agc x = (X 0j i, Y"_i i , M ) C sct) 1 (without time-translations) 
which is a true subalgebra of sci) 1 . Extensions to d > 1 are straightforward. 

What is the usefulness of knowing dynamical symmetries of free, simple diffusion 
for the understanding of non-equilibrium kinetics ? One way of setting up the problem 
would be to write down a stochastic Langevin equation for the order-parameter. The 
simplest case is usually considered to be a dynamics without macroscopic conservation 
laws (model A), where one would have [75] 

where V is the Ginzburg-Landau potential and r] is a gaussian noise which describes 
the coupling to an external heat-bath and the initial distribution of <fi. At first sight, 
there appear to be no non-trivial symmetries, because (fT4"j) cannot be Galilei-invariant, 
because of the noise term 77. To understand this physically, consider a magnet which is a 
rest with respect to a homogeneous heat-bath at temperature T. If the magnet is moved 
with a constant velocity with respect to the heat-bath, the effective temperature will 
now appear to be direction-dependent, and the heat-bath is no longer homogeneous. 
However, this difficulty can be avoided as follows [104] : split the Langevin equation 
into a 'deterministic 7 part with non-trivial symmetries and a 'noise 7 part and then show 
using these symmetries that all averages can be reduced exactly to averages within the 
deterministic, noiseless theory. Technically, one first constructs in the standard fashion 
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(Janssen-de Dominicis procedure) E2] the associated stochastic field-theory with 
action J[<f>, 0] where is the response field associated to the order-parameter 0. Second, 
decompose the action into two parts 

J[0,0] = J o [0,0] + J b [0] (15) 

where 

J o [0,0]=/ dtdr ( 2Md t (j) - A0 + ^ ] (16) 



contains the terms coming from the 'deterministic' part of the Langevin equation (V is 
the self-interacting 'potential') whereas 

J h \$] = -T f dtdr 0(t, rf - \ [ drdr' 0(0, r)a{r - r')0(O, r') (17) 

contains the 'noise'-terms coming from ( 1T41) [82]. It was assumed here that (0(0, r)) = 
and a(r) denotes the initial two-point correlator 

a(r) := C(0, 0; r + r', r') = (0(0, r + r')0(O, r')) = a(-r) (18) 

while the last relation follows from spatial translation-invariance which we shall admit 
throughout. 

It is instructive to consider briefly the case of a free field, where V = 0. Variation 
of (fT5|) with respect to and 0, respectively, then leads to the equations of motion 

2Md t (j) = A0 + T0 , -2Md t 4> = A0 (19) 

The first one of those might be viewed as a Langevin equation if is interpreted as 
a noise. Comparison of the two equations of motion ( fT9j) shows that if the order- 
parameter is characterized by the 'mass' A4 (which by physical convention is positive), 
then the associated response field is characterized by the negative mass —M.. This 
characterization remains valid beyond free fields. 

We now concentrate on actions Jq[0, 0] which are Galilei-invariant. This means 
that if (.)o denotes the averages calculated only with the action Jo, the Bargman 
superselection rules [7] 

~ Sn,m (20) 



hold true. It follows that both response and correlation functions can be exactly 
expressed in terms of averages with respect to the deterministic part alone. For example 
(we suppress for notational simplicity the spatial coordinates) [104] 



0(t)0(,)) = (0(t)0( S )e" J ^) =(0(t)0(s)) (21) 
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\2 ' 

' 



where the 'noise' part of the action was included in the observable and the Bargman 
superselection rule (1201) was used. In other words, the two-time response function does 
not depend explicitly on the 'noise ' at all. The correlation function is reduced similarly 

C{t, s;r)=T [ dudR U(t, r + y)4>(s, y)4>{u, R) 
+ \f dRdR'a(R-R')U(t,r + y)<j>(s,y)0(O,R)0(O,R')) (22) 

l J R 2d \ I 

Only terms which depend explicitly on the 'noise' remain - recall the vanishing of the 
'noiseless' two-point function (0(t)0(s)) o = because of the Bargman superselection 
rule. 

Therefore, the dynamical symmetries of non-equilibrium kinetics are characterized 
by the 'deterministic' part of Langevin equation. Such deterministic non-linear 
diffusion/Schrodinger equations with agc 1 or sci) 1 as a dynamical symmetry can be 
explicitly constructed |119j but we shall not go into the details here. Since all 
quantities of interest will reduce to some kind of response function, one may calculate 
them from the requirement that they transform covariantly under the action ageing 
subgroup (with Lie algebra agc d ) obtained from the Schrodinger group when leaving out 
time-translations. In this survey, we shall concentrate on the two-time autoresponse 
function R(t, s) for which the requirement of covariance reduces to the two conditions 
X R(t, s) = XiR(t, s) = 0. Since time-translations are not included in the ageing group, 
the generators X n can be generalized from (fTTl) to the following form 

X n = -t n+l d t - ^-^t n rd r - (?l±ijl? ' Mt n ~ l r 2 - | (w + l)t n - gnt" ; n>0 (23) 

where £ is a new quantum number associated with the field <fi on which the generators 
X n act. This last term can only be present for systems out of an equilibrium state (the 
requirement of time-translation invariance and [Jfi, Jf_]] = 2X lead to £ = 0). Solving 
the two differential equations for R gives the explicit form of R(t, s), see (126|) below. 

While this discussion was carried out explicitly for the case z — 2, it is tempting 
to try and generalize this idea to more general values of z. In this way, the notion of 
local scale-transformation has been introduced, which is based on the following main 
assumptions [58] . 

(i) In principle, the following conformal time-transformations should be included 

t^t' = ! a8-p>y=\ (24) 

For applications to ageing, however, time-translations generated by f3 must be left 
out (generalizing the restriction sc\) d — > agc d ). 

(ii) The generator X of scale-transformations is 

X = -td t - -rd r - - (25) 

z z 
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Table 1. Magnetic systems quenched into the coexistence phase (T < T c ) which satisfy 
(|26|) with the exponents a — a' and \r. d is the spatial dimension and the numbers 
in brackets estimate the numerical uncertainty in the last digit (s). In the spherical 
model, long-range initial conditions are included and in the long-range spherical model 
the exchange couplings decay as J r ~ |r|~ d ~ CT . In the bond-disordered Ising model, 
the couplings are taken homogeneously from the interval [1 — e/2,1 + e/2}. Then 
z = z{T, e) = 2 + e/T exactly [ElllIOl] and one observes roughly 1.3 < X R (T, e) < 1.7. 



model 


d 


z 


a = a' 






Ref. 


Ising 


2 


2 


1/2 


1.26(1) 




EH 




2 


2 


~ 0.5 


1.24(2) 




[861 EQ] 




3 


2 


1/2 


1.60(2) 




EH 


Potts-3 


2 


2 


0.49 


1.19(3) 




[86, 80J 


Potts-8 


2 


2 


0.51 


1.25(1) 




USED] 


XY 


3 


2 


0.5 


1.7 







XY spin wave 


> 2 


2 


d/2 - 1 


d 


angular response 


[104J 


spherical 


> 2 


2 


d/2 - 1 


(d-a)/2 


Cini(r) ~ k|- d - Q 


[91 [103J 


long-range 


> 2 


cr 


d/a- 1 


d/2 


< a < 2 




spherical 


< 2 


a 


d/<r- 1 


d/2 


< cr < d 




diluted Ising 


2 


2 + e/T 


l/z(T,e) 


\r(T,s) 


disordered 


[72J 



where x is the scaling dimension of the quasi-primary operator on which Xq is 
supposed to act. Physically, this implies that there is a single relevant length scale 
L(t) ~ t 1 / 2 . 

(iii) Spatial translation-invariance is required. 

Generators for infinitesimal local scale-transformations have been explicitly constructed 
and it can be shown that for any value of z there is a linear invariant equation, 
analogous to (jSJ) [58] ■ Local scale-invariance (LSI) assumes in particular that the two- 
time response functions transform covariantly under these local scale-transformations, 
hence X R = X X R = 0. This leads to the prediction [581 (M M, EQ 

R(t, s) = (<Kt)4{ 8 )) = s-^Mt, s) , f R (y) = f y 1+a '- XR/z {y - (26) 

where the exponents a,a',\n/z are related to x, £,x, £ and /o is a normalization 
constant 

Starting with [57], the prediction ff26l) has been reproduced in many different spin 
systems and we list examples quenched to below criticality in table [T] and quenched 

| We point out that the prediction ([26|) as well as the explicit form (|23jl of X n , valid for z — 2, assume 
that the mean order-parameter (0(0, r)) = mo = at the initial moment when the quench to T < T c 
or T — T c is made. 
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Table 2. Systems quenched to a critical point of their stationary state which satisfy 
(I26[) with the exponents a, a' and Xr/z. d is the spatial dimension and the numbers 
in brackets estimate the uncertainty in the last digit (s). CSM stands for the spherical 
model with a conserved order-parameter, FA denotes the Frederikson- Andersen model, 
nekim is the non-equilibrium kinetic Ising model and BCP and BPCP denote the 
bosonic contact and pair-contact processes (see eqs. (|45I48[) for the definitions of 
the control parameter a and of ac), respectively. In the spherical model, long-range 
initial correlations Ci n i(r) ~ |r| -d ~ Q were considered. If d + a > 2, these reduce to 
short-ranged initial correlations (denoted s), but for d + a < 2 a new class L arises. 
In those models described by a Langevin equation, one has used throughout, with the 
exception of the CSM, the simple white noise (n(t, r)i](s, r')) = 2T5(r — r')6(t — s). 



model 


j 
a 


a 


a' — a 






Ref. 


random walk 




-1 










[34j 


OJK-model 




(rf-l)/2 


1 /O 

-1/2 




a/4 


PZIGIHE21 


Ising 


-i 

1 





1 /o 

-1/2 




1 /o 

1/2 




[501 [851 [62] 




o 
I 


n iir 


n i 07/ on\ 
— U.lo^zUJ 




n 'TOO / cr \ 

U. Mz^oJ 




|108[ E] 




Q 
O 


n ^dr 






l.OU 




[1081171] 


XY 


3 


0.52 







1.34(5) 





spherical d > 2 


< 4 


d/2 - 1 





d/2- 


-a/4-1/2 


L 


[103] 




> 4 


d/2 - 1 





(d- 


a)/4+ 1/2 


L 


[103] 




< 4 


d/2 - 1 







3d/A - 1 


s 


EU 




> 4 


d/2 - 1 







d/2 


s 


[51] 


CSM 


> 2 


d/A-1 







(d + 2)/4 


[14j 


disordered Ising 


4-e 


1 1 h e 

2 y 53 







1 /6e 
° 2 V 53 


0(e), log 


|23L 11151 1116] 


FA 


> 2 


l + d/2 


-2 




2 + d/2 




m 




1 


1 


-3/2 




2 




[9Ql[88] 


Ising spin glass 


3 


0.060(4) 


-0.76(3) 




0.38(2) 


[66, 69J 


contact process 


1 


-0.681 


+0.270(10) 




1.76(5) 


t/s> 1.1 


[HI [751 nu 




> 4 


d/2 - 1 







d/2 + 2 




[inn] 


NEKIM 


1 


-0.430(4) 







1.9(1) 


[96j 


BCP 


> 1 


d/2 - 1 







d/2 


IS] 


BPCP 


> 2 


d/2 - 1 







d/2 


a < «c 






to the critical point in table [3 For T < T c , it is found empirically that a = a' in all 
examples considered so far. We point out that agreement with local scale-invariance 
eq. (|26l) is not only obtained for systems where the dynamical exponent is z = 2, but 
that rather there exist quite a few examples where z can become considerably larger or 
smaller than 2. It must be remembered, however, that the above derivation of f[26l) for 
a stochastic Langevin equation has for the time being only been carried out for z = 2^| 



§ See section 4 for a recent extension of the method to z = 4. 
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and the justification of X R = XiR = remains an open problem for z 7^ 2 although 
the result fl26l) seems to work remarkably well. Still, it is non-trivial that a relatively 
simple extension of dynamical scaling should be capable of making predictions which 
can be reproduced in physically quite different systems. 

A few comments are still needed: (i) for the XY model in the spin-wave 
approximation (tabled]), eq. (1261) holds for the response of the angular variable <fi = 
4>(t,r) which is related to the XY spin through S = (cos <ft, sin <ft). Magnetic responses 
have a different scaling form [18j[l]. (ii) in the critical disordered Ising model (table W) 
one finds a logarithmic scaling form R(t, s) = (r + 7*1 ln(£ — s))/i?(i/s) [231 IH5L 1116] 
such that the computed fniy) is consistent with (1261) to one- loop order, or up to terms 
of order 0(e). (iii) Finally, a two-loop calculation of the critical non-conserved O (re- 
model does produce in 4 — e dimensions an expression for fn(y) which is incompatible 
with ( 1261) [22] and a similar result is anticipated in 2 + e dimensions [15]; although 
the one-loop results are still compatible [211 1221 H5]- Should one conclude from these 
studies that for T = T c the prediction (|26|) and by implication local scale-invariance 
can only hold approximatively ? This might well be a subtle question. Deviations 
between fl26l) and the field-theoretical studies typically arise when t/s ~ 1. However, 
in this region the field-theoretical results for /#(?/) do not agree with the ones of non- 
perturbative numerical studies |108| . Since the perturbative expansion usually carried 
out in field-theoretical studies does not necessarily take care of the Galilei-invariance, it 
is necessary to carefully check that the truncation of the e-series does not introduce slight 
inaccuracies. Only after this has been done (for example by re-summing the e-series) and 
checked by comparing with non-perturbative data, meaningful quantitative statements 
on the scaling functions can be made, (iv) Throughout, it was implicitly assumed that 
the order-parameter vanishes initially. Systematic studies on what happens when this 
condition is relaxed are only now becoming available [51 [251 H5] • These extensions might 
be particularly important for chemical kinetics, (v) We did not include growth models 
here but shall discuss them in section 4. 

If z = 2, it is also possible, using eq. ( 1221) . to derive explicit predictions for the 
two-time correlation function [1041 [67] . These have been tested in some exactly solvable 
models [104, 71J, the 2D Ising model [67] and the 2D g-states Potts model with q = 2, 3, 8 
|86[ [80] . Extensions to z = 4 have been studied very recently [1101 fT4"] - see section 4. 

This survey is organized as follows. In section 2 we review results on the ageing 
behaviour of several critical models with non-equilibrium steady-states. The first two 
models are chosen because their steady-state phase-transitions are in the paradigmatic 
universality classes of directed percolation (DP) and in the parity-conserving (PC) 
universality classes. The numerical results on these models are supplemented by the 
exactly solved bosonic variants of the contact and pair-contact processes. In this way it 
becomes clear that most aspects of the scaling description of ageing in simple magnets 
does carry over to this more general situation. However, a central issue, namely the 
definition of an universal limit fluctuation-dissipation ratio X^ [51] and which has 
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received so much recent attention in the non-equilibrium critical dynamics of magnets, 
see j32j [24] for reviews, does not appear to have an obvious analogue. Remarkably, 
the same kind of evidence in favour of a non-trivial extension of dynamical scaling 
towards a larger dynamical symmetry group of local scale-transformation previously 
found in magnets also appears in the models without detailed balance. In section 3, we 
review in more detail how the stochastic Langevin equations underlying the bosonic 
contact and pair-contact processes can be shown to actually possess a local scale- 
invariance. These evidences should form a promising basis to look for more manifestation 
of local scale-invariance in systems with dynamical scaling which remain always very 
far from equilibrium. A different class of non-equilibrium models is studied in section 4, 
where kinetic growth as described by the Edwards- Wilkinson and the Mullins-Herring 
equations is studied. To account for those models, the formulation of LSI must be 
generalized [58J to values z ^ 2 of the dynamical exponent. We consider explicitly the 
case z = 4 and apply it to the Mullins-Herring model. We conclude in section 5. 

2. Ageing with absorbing steady-states 

We now describe the ageing behaviour of system without an equilibrium stationary 
state. We shall realize this system as models of interacting classical particles, where the 
stochastic dynamics is such that the detailed-balance condition no longer holds. For the 
simple models we shall consider here it turns out that if the stationary state is not at a 
critical point, only a single stable stationary state remains to which the system relaxes 
within a finite time and no ageing is possible. For this reason, we shall study the ageing 
behaviour at criticality. 

We remark that one might also go to non-equilibrium stationary states by 
considering driven systems [1 13j . However, the dynamics of those is more complicated 
than the systems at hand because of a further strong spatial anisotropy and the 
description in terms of local scale-invariance would require to generalize the local scale- 
transformations accordingly. That is beyond the scope of this survey. Another very 
interesting class of ageing non-equilibrium systems are zero-range process, see [431 [53] 
for recent reviews. Because they do not have a spatial structure which would admit 
a Galilei-invariance, their dynamical scaling cannot be extended to some form of local 
scale invariance and for this reason they are not considered here, despite their intrinsic 
interest. 

2.1. Contact process 

The contact process is a paradigmatic system for the study of non-equilibrium phase- 
transitions, see e.g. [74] for a review [jj] The steady-state phase-transition of the contact 

|| Recall that bifurcations arising in many simple models of mathematical biology [93j come from the 
mean-field treatment of the phase-transition in the contact process. 
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Figure 3. Microscopic evolution of clusters in the critical 2D contact process, on a 
lattice of size 1000 x 1000. The initial condition of the upper series is a full circle with 
radius 100 placed in the center of the lattice, while in the lower series it is a full lattice. 
The times are t = [2, 20, 200, 2000] for the upper series and t = [20, 200, 2000, 20000] 
below. After [109]. 



process is in the same universality class as one of the transitions of the celebrated 
Ziff-Gulari-Barshad model [128] . which is meant to describe the catalytic reaction 
2CO + O2 — > 2CO2. The model may be defined in terms of a time- dependent discrete 
variable rii(t) G {0, 1}, defined on each site i of a hypercubic lattice, which describe 
configurations of particles and empty sites. The dynamics is defined as follows: for each 
time-step, select randomly a site i of the lattice. If i is occupied (i.e. = 1), that 
particle vanishes with probability p. Otherwise, with probability 1-pa new particle 
is created on one of the nearest neighbours of i, chosen at random and provided that 
chosen site is still empty. Formally, this may be expressed through the reactions A — > 
and A — ► 2 A, with rates corresponding to p and 1— p, respectively. In the steady-state, 
the model has a continuous phase-transition at some critical value p c . Numerically, 
p c = 0.2326746(5) in ID and p c = 0.37753(1) in 2D. 

A first characteristic of the dynamics of the critical contact process can be seen 
by looking at the temporal evolution of certain initial configurations, see figure [3j 
In contrast to magnetic systems, see figure [2] for comparison, in the contact process 
there is no apparent growing length-scale at all and the evolution proceeds via the slow 
dissolution of the particle clusters. Cluster dilution had first been demonstrated to occur 
in several variants of the two-dimensional voter model [38] but also occurs in the early 
stages of surface ageing in simple magnets [1061 1107] . 

In studies of the ageing behaviour, one goes beyond the average particle-density 
N(t) := (ni(t)) ~ t~ s at criticality p = p c . Define the two-time (connected) 
autocorrelator and autoresponse functions 



C(t,a):= <m(f)"i(a)> - <n,(t)>Ma)) , R(t,s):- ' H "' (/)> 



Shi(s) 



(27) 



h=o 

where hi(s) is the rate of the spontaneous creation process — > A at the site i at time 



Ageing & scaling without detailed balance 



15 




Figure 4. Connected autocorrelation function C(t,s) and autoresponse function 
R(t,s) of the ID contact process in the active phase (p = 0.1). The straight lines 
are proportional to exp(— 0.05(f — s)). After [41] . 

s. These may be calculated in a standard fashion either from simulations [109] [75] 
or else from the transfer-matrix renormalization group [41] . Recent field-theoretical 
calculations p3] will also be described. 

2.1.1. Active phase In contrast with simple magnets, where there are two distinct 
stable ground states in the low-temperature phase, in the active phase of the contact 
process there is only a single stable steady state. Consequently, there is here no breaking 
of time-translation invariance and we illustrate this in ID in figure HJ After a short 
transient, the data for both C(t, s) and R(t, s) collapse when plotted over against t — s 
which means that the contact process shows no ageing in its active phase. 

2.1.2. Absorbing phase From the comparison with the high-temperature phase of 
simple magnets, one would also expect to find time-translation invariance in the 
absorbing phase of the contact process. However, the correlation function shows a 
subtlety the origin of which is best understood by considering the case p = 1 first. If 
p = 1, particles on different sites are uncorrelated and simply decay with a fixed rate. 
For any fixed site % and with two times t > s, it is clear that ni(t)rii(s) = n^t), since 
ni e {0,1}. Hence (m(t)m(s)) = (rij(t)) and C(t,s) = N(t)(l - N(s)). For sufficiently 
long times, C (t, s) will then only depend on t. Indeed, this behaviour survives in the 
entire absorbing phase [41J. On the other hand, the expected time-translation invariance 
for the autoresponse function is readily checked. 

2.1.3. Critical point For the critical contact process we show in figure [5|l that ageing 
does occur, that is, the autocorrelation and the autoresponse depend on both the 
observation time t and the waiting time s. Furthermore, when the same data are 
re-plotted over against t/s, a data collapse after rescaling can be achieved, see figure [5b. 
Lattices with a large initial particle density n « 0.8 — 1 were used [41], 11091 EH]. This is 
different with respect to the magnetic systems of section 1, where the order-parameter 
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Figure 5. Connected autocorrelation function of the critical contact process in ID 
(main plots) and 2D (insets). Panel (a) shows the ageing of the autocorrelation 
functi it lines correspond 

to the I I 




t/s 



Figure 6. Autoresponse function for the critical ID contact process for several waiting 
times s. The data labelled tm come from the transfer matrix renormalization group 
Pffl and MC denotes Monte Carlo data. 



had a vanishing initial value. By analogy with simple magnets, one defines the ageing 
exponents a, b and the autocorrelation and autoresponse exponents Xc,r from 

C(t,s) = s- b f c (t/s) , f c (y)~y~ Xc/z 

R(t, s) = s-^Mt/s) , f R (y) ~ (28) 

where the asymptotic forms should hold for y — ► oo. Similarly, scaling can be observed 
for the autoresponse function as shown in figure El On the other hand, the unconnected 
correlator behaves for large times simply as (ni(t)n,i(s)) ~ (ts)~ b / 2 . 

The results for the ageing exponents a, b, Ac, Xr are collected in table [3j The 
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Table 3. Nonequilibrium exponents for the contact process (cp), the non-equilibrium 
kinetic Ising model (nekim), the bosonic contact process (bcp) and the bosonic 
pair-contact process (bpcp). Several kinetic growth models based on the Edwards- 
Wilkinson (ew) and Mullins-Herring (mh) equations are also listed, see section 4. 
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agreement between the results of the TMRG and Monte Carlo simulations serves as 
a useful controle on the fiability of the results. While the equality Ac = \r is fully 
analogous to what was seen in non-equilibrium critical dynamics, the exponents a and 
b are no longer equal but satisfy 

1 + a = b = 25 (29) 

Some comments are in order. 



(i) If the critical contact process were a Markov process, these exponents might be 
calculated from the global persistence exponent 9 g through the scaling relation 

(E7J [73113] 

hi = tt _ HClz^ (30) 

which would predict Xc/z = 1.98(2) in ID and Xc/z = 3.5(5) in 2D. Although 
this not too far form the values reported in table [31 the differences appear to be 
significant. If that conclusion is correct, it would point towards the existence of 
temporal long-range correlations and hence of an effective non-markovian dynamics 
of the critical contact process. 

ii) The relation 1 + a = b can be understood [13] as follows to be a consequence of 
the rapidity-reversal symmetry of Reggeon field-theory (which is generally thought 
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to be in the same universality class as the contact process). The field-theoretical 
action in the Janssen-de Dominies formulation reads at the critical point 

J[<j>,j>] = Jdtdr (fi(d t — DA) <j) — u (f?J (f)(f) — h4> (31) 

where <fi and h are the coarse-grained particle-densities and creation rates for 
particles. For h = and if the time i 6 K is unbounded, the action is invariant 
under the rapidity reversal, see |121l I122j 

0(t, r) ^ -0(-t, r) , d)(t, r) ^ -0(-t, r) (32) 

In particular, it follows that the scaling dimensions x^ = x^ = f3/v must be equal. 
This remains true even if rapidity-reversal is broken by initial conditions at time 
t = 0. For a rapidity-reversal-invariant action J the connected correlator becomes 



C(t, s;r- r') = (<Kt, r)0(s, r')> - (<f>(t, r))(<f>(s, r')) 

= (<j>(-t, r)<K-s, r')> - (0(-t, r)) (</>(- s, r')) 

= (33) 

where the second line comes from rapidity-reversal symmetry and the last line 
follows from causality. Hence C(t, s; r) = in the steady-state but for relaxations 
from an initial state C(t,s;0) = (4>(t)4>(s)) c and R(t, s; 0) = (<p(t)4>(s)) are non- 
vanishing and have the same scaling dimensions, which implies (J29l) and also 
Ac = ^R- 

While this explains the origin of eq. (|29|) for the contact process, it is not yet 
understood why it also holds true in the nekim [96] (see below) where rapidity- 
reversal symmetry is not known to be satisfied. On the other hand eq. (1291) is not 
universally valid. In the critical bosonic contact process, one has a = b. The critical 
bosonic pair-contact process furnishes further examples with a ^ b but the relation 
between a and b is distinct from eq. (1291) . 

(iii) The non-equality of the exponents a and b is only possible in systems with non- 
equilibrium steady-states. Indeed, for equilibrium systems, one has time-translation 
invariance and combining this with the scaling forms would give 

C{t,s) ~ {t-s)~ b , R{t,s) ~ {t-s)- 1 -' 1 

The fluctuation-dissipation theorem would then give a = b. Hence the equality 
a = b is a necessary condition that a quasi-stationary state, which might be present 
for t — s <C 1, is an equilibrium one. 

(iv) Is it possible to define a non-equilibrium temperature from the steady-states of 
systems without detailed balance ? A recent attempt by Sastre, Dornic and 
Chate [112] started from the observation that in simple magnets, the fluctuation- 
dissipation ratio X(t, s) — > 1 as t — > oo and t — s — > 0. From this observation, they 
define a dynamical temperature by 

1 ■ lim ( hm ) (34) 



T dyn ' t^oo \t-s-* dC(t, s) /ds 
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By explicit calculation, they confirm that in the 2D critical voter model (where 
indeed a = b) this limit exists, has a non-trivial value and is universal [1 12j . Still, 
their appealing idea has met with several criticisms. First, Mayer and Sollich 
[89] construct in the coarsening ID Glauber-Ising model a defect-pair observable 
such that the fluctuation-dissipation ratio X(t, s) ^ 1 in the short time- regime (in 
particular they show lim^oo X(s, s) = 3/4). Second, T dyn can only be finite if a = b 
and the examples listed in table [3] show that the fluctuation-dissipation ratio itself 
is in general no longer defined. It appears that the proposal ( 1341) relies too heavily 
on specific properties of the voter model. 

Rather than the fluctuation-dissipation ratio X(t, s) as defined in eq. ([3]) for systems 
with detailed balance, one may instead consider the ratio E(t, s), and its limit H^, 
which are defined by [JT] 

R(t, s; 0) f R (t/s) „ ( \ 

which are well-defined because of ( 1291) . The limit 2^, being the ratio of two 
quantities with the same classical and scaling dimensions, is expected to be universal 
[T3] . Baumann and Gambassi [13] argue that a value E(t, s)^ 1 ^ is a measure 
for the breaking of the rapidity-reversal symmetry (in the same way as X ^ 1 
measures the distance from an equilibrium state) and show explicitly that only the 
zero-momentum modes contribute to a non- vanishing value of S(t, s) _1 . In this 
respect, for systems with rapidity-reversal symmetry, S appears to be the analogue 
of the fluctuation-dissipation ratio R. Of course, the applicability of their argument 
depends on the validity of the scaling relation 1 + a = b. For the value of Soo, they 
find in 4 — e dimensions 

119 7T 2 



1 - e — - — + 0(e 2 ) (36) 

V480 120 y J v ; v ; 

which in ID is in semiquantitative agreement with the estimate Soo = 1.15(5) [41] . 

Before one can discuss the form of the scaling function fn(y), it is necessary to study 
the role of the initial conditions. Indeed, all existing simulations on the ageing in that 
model start from a lattice with a non-vanishing particle density and hence in contrast 
to simulations in magnetic system, where the initial order-parameter was set to zero. 
For the time-dependent order-parameter one expects the scaling form [1261 1 



W)) = 4>ot e f {U e+P,(vz) ) (37) 

where 0o = (0(0)) is the initial value of the order-parameter (which for the contact 
process is the particle density) and 9 is the slip exponent. Hence there is characteristic 
time scale r* ~ <p 1 ^ e+l3 ^ uz ^ where a change of scaling behaviour takes place. For any 
non-vanishing value of the dimensionful variable 0o the long-time scaling behaviour is 
effectively described by the r* — > limit. For a vanishing initial order-parameter the 
exponents Xc,r are independent of the equilibrium exponents [ST] . For O 7^ 0, Baumann 
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and Gambassi show from an analysis of the scaling behaviour of both responses and 
correlators that [13] 

\ c = \ R = d + z + ^ (38) 

The available exponent estimates in ID and 2D from simulational studies [ITj 11091 I75| 
ITT] agree quite well with this prediction. We point out that the available results for the 
limit ratio [UJ |J3], see eq. are found for a finite initial particle-density. 

On the other hand, in the theory of LSI as reviewed in section 1 the implicit 
assumption <p = was made. Since O couples to a relevant scaling variable, it is not 
obvious why a theory formulated in the limit (p t e+l3 ^ uz ^ — > should be valid in the 
opposite limit <j) t e+l3 ^ uz ^ — > oo where all existing studies on the ageing in the contact 
process have been performed. It is therefore remarkable that LSI with 0o — still 
captures well the behaviour of the 0o 7^ data. This is shown in figure El where 
for almost all values of t/s an almost perfect agreement with the LSI prediction (j26|) . 
derived for O = 0, is foundQ Hinrichsen [75] has argued that a more ambitious test 
may be performed by plotting r(y) := fR(y)[y XR ^ z ~ 1 ~ a ' (y — l) 1+a '] over against t/s — 1. 
If LSI with 0o — holds, one expects that r(y) = const. In this way, he found that for 
although the measured function fn(y) agrees nicely with eq. (126]) if t/s is large enough, 
for values t/s < 1.1, the function fii(y) remains well-defined but changes to a different 
behaviour which is no longer described by eq. ( 1261) . The high quality of his data makes 
it clear that this change of behaviour in the scaling function cannot be explained away 
by invoking corrections to scaling. Further unpublished calculations [12] for extremely 
large values of s confirm these conclusions. In addition, the same conclusion has been 
reached from a detailed field-theoretic study of the two-time response in momentum 
space [13]. A possible explanation of the form of /#(?/) in terms of LSI will require the 
extension of the theory to include a non- vanishing value of 4>q. 

2.2. Non- equilibrium kinetic Ising model 

We now review results, obtained by Odor |96j, on the ageing behaviour in a non- 
equilibrium kinetic Ising model (nekim) where the parity of the total particle-number 
is conserved. The model was introduced by Menyhard [92] and combines spin-flips as in 
zero-temperature Glauber dynamics with spin-exchanges as in Kawasaki dynamics. The 
model is formulated either in terms of Ising spins (f, |) or else in terms of a particle- 
reaction model of the kinks with occupied or empty sites (•, o). First, the Glauber-like 
part of the dynamics contains a diffusive motion 

TII^TTI or equivalently • o ^ o« ; with rate D 
In the limit <f>o — > 0, the results for R(t, s) appear to be consistent with LSI [13] . 
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and the pair-annihilation of nearest neighbours 

tit— >ttt or equivalently • • — > oo ; with rate 2a 

The Kawasaki-like part of the dynamics is described by 

TTII^TITI or equivalently o mo ^± m m m ; with rate k 

In full, this is a model describing branching and annihilating random walks with an 
even number of offspring. By increasing k, one finds a second-order phase-transition 
[92] where the kinks go from an absorbing to an active state. This phase-transition is 
in the so-called parity-conserving (PC) universality class [5"H [29| [26] which is different 
from the one of the contact process. 

Using the parameterization k = 1 — 2T, D = r(l — S)/2 and 2a = T(l + S), the 
critical point is located at T — 0.35, k = 0.3 and 5 = —0.3928 [96]. Measuring the 
kink-density through an efficient cluster algorithm, and starting from a fully ordered 
kink state with spins being alternatingly j and j, he finds a nice power-law scaling 

(m(t)} ~r°- 285 ( 2 ). 

Next, measuring the unconnected kink-kink two-time correlation function, Odor's 
data are fully compatible with the scaling behaviour (ni(t)ni(t')) ~ t' _0 - 57 (t/t / ) -0 - 285 for 
t/t' — > oo and allows to determine the ageing exponent b, see table [31 This result is 
also consistent with earlier results on the spin-spin autocorrelator in the same model, 
see |96j for details. The connected autocorrelator was also calculated, with the result 
Xc/z = Xr/z = 1.9(1). This is in agreement with the scaling relation (|38|) . 

Finally, the spin-autoresponse R(t, s) with respect to a magnetic field coupling to 
a spin was calculated [96], by adapting methods used previously by Hinrichsen [73] for 
the contact process. We recall that for the fully ordered kink initial state, the initial 
spin magnetization vanishes. In figure [7j it is shown that a clear data collapse is found 
when 1 + a = b, as in the contact process. 

Comparing the form of the scaling function fB,{y) with the prediction of LSI, Odor 
finds a perfect agreement as long as his data remain in the scaling limit. In particular, 
he carefully considered the limit y — > 1 (see inset in figure [7J). Within the numerical 
accuracy, the data for the rescaled response scaling function r(y) remain essentially 
constant for all values oit/s, as long the model is in the scaling regime. Clear deviations 
from a horizontal line only occur once finite-time corrections have broken dynamical 
scaling and by increasing s, the regime where (i) scaling holds and (ii) r(y) is constant 
progressively extends to ever smaller values of y = t/s. This finding is different from 
the ID contact process [751 EE] ( see above). It is conceivable that the agreement with 
LSI in the nekim comes from the fact that the initial configuration used has a vanishing 
magnetization. 
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2.3. Bosonic contact and pair-contact processes 

Two exactly solvable models allow to confirm the conclusions drawn from the above 
numerical studies. 

The bosonic contact process was introduced in order to describe clustering 
phenomena in biological systems [77] whereas the bosonic pair-contact process was 
originally conceived [HS] as a solvable variant of the usual (fermionic) pair-contact 
process. These models are defined as follows. Consider a set of particles of a single 
species A which move on the sites of a rf-dimensional hypercubic lattice. On any site one 
may have an arbitrary (non-negative) number of particles and it is this property which 
makes up the difference with the usual contact and pair-contact processes considered 
before where on each site at most one particle is allowed. Single particles may hop to a 
nearest-neighbour site with unit rate and in addition, the following single-site creation 
and annihilation processes are admitted 

mA — (m + 1)A , pA — (p — t)A ; with rates \i and A (39) 

where t is a positive integer such that \l\ < p. We are interested in the following special 

cases: 

(i) critical bosonic contact process: p — m — 1. Here only i = 1 is possible. 
Furthermore the creation and annihilation rates are set equal fi = A. 

(ii) critical bosonic pair-contact process: p = m = 2. We fix I = 2, set 2A = \x and 
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define the control parameter 

* := % (40, 

The dynamics of these models is conveniently described in terms creation operators 
a){t, r) of a particle at time t and location r and the corresponding annihilation operator 
a(t,r), see [36J 1117] . The equation of motion for the space-time-dependent particle- 
density p(t, x) := (a)(t, x)a(t, x)) = (a(t, x)) reads, after a rescaling t \—> t/(2D) [771199] 

^ (a(t, x)) = U x (a(t, x)) - ^ (a(t, xf) + JL (a(t, x) m ) + h(t, x) (41) 
where we have used the short-hand 

d 

A x f(t, x) := (fit, x-e r )+f(t,x + e r ) - 2f(t, x)) (42) 

r=l 

Similarly, the equal-time correlation functions satisfy the equations of motion 



d_ 

di 



1 d 

(a(x)a(y)) = - ^ \(a(x)a(y - k)) + (a(x)a(y + k)} 



^ ((a(x)) 2 ) = [«x)a(x - fe)> + (o(aj)o(a! + fc)> - 2 (a^) 2 ) 

k=l 



k=l 

+ (a(x - fc)a(y)) + (a(» + fc)a(y)> - 4 (a(aj)a(y)) ( 13) 

d 



+ HxD (44) 

where in eq. (j43p x ^ y is understood. Since (n(x) 2 ) = (a(x) 2 ) + (a(a;)), the main 
equal-time quantity of interest, namely the variance a 2 := (n(x) 2 ) — (n(x)) 2 can be 
found. 

The equations of motion f)43|44l) are already written for the critical line given by 



t\ = p. (45) 

where they naturally close. For the bosonic contact process p = m = I there is an 
extension to arbitrary values of A, /i which still closes. In both models, the spatial 
average of the local particle-density p(x,t) := (a(x,t)). remains constant in time 

Jdxp(x,t) = J dx (a(x,t)) = p (46) 

where po is the initial mean particle-density. Furthermore, the critical line (l4"5j) separates 
an active phase with a formally infinite particle-density in the steady-state from an 

+ If instead we would treat a coagulation process 2 A — > A, where i = 1, the results presented in the 
text are recovered by setting A — fi and a = 
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absorbing phase where the steady-state particle-density vanishes. The phase diagrams 
are sketched in figure 



X 




X 




(b) 



Figure 8. Schematic phase-diagrams for D ^ of (a) the bosonic contact process and 
the bosonic pair-contact process in d < 2 dimensions and (b) the bosonic pair-contact 
process in d > 2 dimensions. The absorbing region 1, where limt_ >oo p(x, t) = 0, is 
separated by the critical line eq. from the active region 2, where p(x,t) — > oo as 
t — ► oo. Clustering along the critical line is indicated in (a) and (b) by full lines, but in 
the bosonic pair-contact process with d > 2 the steady-state may also be homogeneous 
(broken line in (b)). These two regimes are separated by a multicritical point. 



The physical nature of this transition becomes apparent when equal-time 
correlations are studied [77J ES]- For example, for the bosonic contact process at 
criticality one has [77] 



(a(t,x) 2 ) 



a t~ d / 2+1 

C2 hit 

c 3 + c 4 t~ d / 2+1 



if d < 2 
if d = 2 
if d > 2 



(47) 



where f > 1 and cq, . . . , C4 are known positive constants. For d < 2, the fluctuations in 
the mean particle-density increases with time, although the mean particle-density itself 
remains constant. Physically, this means that the particle- number on relatively few sites 
will increase while many other sites will become empty. Only for d > 2 fluctuations will 
eventually die out. For the bosonic contact process, this critical behaviour is the same 
along the entire critical line. 

For the bosonic pair-contact process, that is different. Rather, there exists a critical 
value «c of the control parameter, given by 



1 

a c 



POO 

2 du (e- Au I Q {Au)y 
Jo 



(4£ 



and where Iq(u) is a modified Bessel function. Specific values are «c(3) ~ 3.99 and 
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a c{4) ~ 6.45 and lim(j\2 occ(d) = 0. Then the variance behaves as 



(a(t,x) 2 ) 



fo 

ht 



if a < ac 

if a = ac and 2 < d < 4 
if a = ac and d > 4 
if a > ac (or d < 2) 



(49) 



where the fo, . . . , f% and r are known positive constants. This means that at the 
multicritical point at a = ac there occurs a clustering transition such that for a < « c 
the systems evolves towards a more or less homogeneous state while for a > ac particles 
accumulate on very few lattice sites while the other ones remain empty. In contrast with 
the bosonic contact process, clustering occurs in some region of the parameter space for 
all values of d. 

We are interested in studying the impact of this clustering transition on the two- 
time correlations and linear responses. In order to obtain the equations of motion of 
the two-time correlator, the time-ordering of the operators a(t, x) must be taken in 
account. This leads to the following equations of motion for the two-time correlator, 
after rescaling the times t h- > t/(2D), s \— > s/(2D), and for t > s, [491 15] 

^(a(t,x)a(s,y)) (50) 
= (<*(t, x)a{s, y)) - ^ (a(t, x) p a(s, y)) + ^ (a{t, x) m a{s, y)) 



We are interested in the two-time connected correlation functio 



rS 



C(t, s; r) := (a(t, x)a(s, x + r)) 



(51) 



and take an uncorrelated initial state, hence C(0, 0; r) = 0. The linear two-time response 
function is found by adding a particle-creation term Y2 x h(x,t) [oj\x,t) — l) to the 
quantum hamiltonian H and taking the functional derivative 



R{t,s;r) 



8{a(t, r + x)) 



5h(s, x) 



(52) 



h=0 



and for which the usual scaling behaviour ( 1251) is anticipated. 

The solution of the equations ( 1501) for the two-time quantities is now 
straightforward, if just a little tedious. It can be shown [8] that the anticipated scaling 
behaviour (1251) exists along the critical line, but for the bpcp the further condition 
a < ac is required. For these cases, the exponents are listed in table [3j In particular, 
we see that at the multicritical point a = ac, the exponents a and b are different and 
furthermore, do not satisfy the relation 1+a = b found for the critical contact process and 
the nekim. This means that there is no straightforward way to define an analogue of a 

* It can be shown that in the scaling regime C(t, s; r) ~ r )n(s, r + r )) describes the two-time 
density-density correlation function |12j . 
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fn(y) 


/o(y) 


contact process bcp 


(v-i)-* 




pair 
contact 
process 


a < etc 


d>2 


(v-i)"* 


(„ _ (y + 


a — ac 


2<d<4 


{y-i)-i 




d>4 


(y-i)-5 


(y + l)-i+ 2 - (y - l)-i+ 2 + (d - 4)(y - l)-§ +1 


Edwards- Wilkinson ew2 


(v-i)-* 


(y_ l)-f _ (y + l)-|+l+P 


Mullins-Herring MHl 


(v-i)"* 


(y _!)-!+! _(y + !)-!+! 


Mullins- Herring MH2 


(v-i)-* 


(y _ l)-|+l+p/2 _ ( y + !)-| + l+p/2 


Mullins-Herring MHc 


(y _ l)-(<i+2)/4 


(y _ l)-(^-2)/4 _ (y + ^-(d-^ 



Table 4. Scaling functions (up to normalization) of the autoresponse and 
autocorrelation of the critical bosonic contact and bosonic pair-contact processes [5] . 
The scaling functions for some simple growth models [1101 [T4] . as defined in section 4, 
are also listed. 



limit fluctuation-dissipation ratio for particle-reaction models without detailed balance. 
Furthermore, the explicit form of the scaling functions can also be found and are listed in 
table HI While the form of the autoresponse function fn(y) — (y — l)~ d / 2 is remarkably 
simple, the results for the autocorrelation function can be rendered as an integral 

fc(y) = Co [ d6 6 a -\y + 1 - 2dy d ' 2 (53) 
Jo 

where the exponents a, b are taken from table [31 

In section 3, we shall show how these results for fn(y) and fc(y) in the bcp and 
the bpcp can be understood using local scale-invariance. In section 4, we shall define 
the EW and MH models and, after having briefly reviewed LSI for z — 4, then perform 
a similar analysis. 

3. The bosonic processes and local scale-invariance 

We now show that the exact results for response and correlation functions of the bcp 
and the bpcp as listed in table H] can be understood from local scale-invariance [9l [H] . 

3.1. Bosonic contact process 

The master equation which describes the critical bosonic contact process can be turned 
into a field-theory in a standard fashion through an operator formalism which uses a 
particle annihilation operator a(t,r) and its conjugate a^(t,r) [361 1121 j . For calculating 
connected correlators, it is useful to define the shifted fields 

(f>(t,r) := a{t,r) - p 
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field 


scaling dimension 


mass 






X 


M 







X 


-M 






x 2 


-2M 


T 


:=0 2 


it 


-M 


S 


:= 3 


Xt, 


-2M 


r 


= 3 2 


x r 


-M 



Table 5. Scaling dimensions and masses of some composite fields. 



4>(t, r) := a(t, r) = a\t, r) - 1 (54) 

such that (<f>(t, r)) = (our notation implies a mapping between operators and quantum 
fields, using the known equivalence between the operator formalism and the path-integral 
formulation [521 IT5| ll^lj ) ■ As we shall see, these fields <ft and <p W1U become the natural 
quasi-primary fields from the point of view of local scale-invar iance. We remark that 
the response function is not affected by this shift, since 



R(t,s;r,r') 



6{a(t,r)) 



5h(s,r') 



Sm,r)) 



h=0 



5h(s,r>) 



h=0 



As for magnets, the field-theoretical action |78j is again decomposed J[ 
Jb[4>) 0] into a 'deterministic' part 



Jn 



Jo 



dR du 



(t>{2Md u - V 2 



and which is manifestly Galilei-invariant, whereas the 'noise' is described by 



li dR du 2 (0 + p o ) 



Jn 



We use uncorrelated initial conditions C(0, 0; r) = throughout. 



(55) 



+ 



(56) 



(57) 



In what follows, some composite fields will be needed, which we list, together with 
their scaling dimensions and their masses, in table [5j We remark that for free fields one 
has 



x 2 = 2x , xy = 2x + x , is = 3x + x , xr = 3x + 2x 



(58) 



but these relations need no longer hold for interacting fields. On the other hand, from 
the Bargman superselection rules we expect that the masses of the composite fields as 
given in table [5] should remain valid for interacting fields as well. 

As in section 1, we now have a similar reduction to averages of the noiseless 
theory. First, for the computation of the response function, we add the term 
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J dR J du<p(u, R)h(u, R) to the action. As usual the response function is [9J 
R(t,s ] r,r') = ^(t,r)^(s,r / )^ 

= (${t, r)0(s, r') exp ^-^ J dR j du 4> 2 {u, R)(<p(u, R) + p )^ 

= (^(t, r)0(s, r')) o = i2o(t, s; r, r') (59) 

where we expanded the exponential and applied the Bargman superselection rule. 
Indeed, the two-time response is just given by the response of the (gaussian) noise- 
less theory. We have therefore reproduced the exact result of table 0] for the response 
function of the critical bosonic contact process. 

Second, we have for the correlator 

C(t,s;r,r r ) = ($(t, r)<f)(s, r') exp j dR j 'du4> 2 (u, R)<p(u, R) 



x exp ^— fipo J dR J du (f) 2 (R,u] 



(60) 



o 



Expanding both exponentials separately and using the Bargman superselection rule (1201) 
it follows that C = C\ + Ci can be written as the sums of two terms which read 

Ci(t, s; r, r') = - m J dR J du r )<K s ' r ')4> 2 (^ #)} Q (61) 



and 

,2 



C 2 (t, s; r, r') = Y J dRdR J dudu' (<f>(t, r)<j>(s, r')T{u, R)T{u', R')) (62) 

using the field T, see table Hence the connected correlator is determined by three- 
and four-point functions of the noiseless theory. 

The noiseless three-point response needed for C\ can be found from its covariance 
under the ageing algebra [9] 

(0(t, r)(j)(s, r')4> 2 {u, R)) = (t - s) x ~& (t - u)-^ 2 (s - u)~^ 
Mir- R? M (r' - R^ 2 



with 



x exp ~ — V 3 (u u v l )Q(t - u)<d(s - u) (63) 

2 t — u 2 s — u 

u [(s — u)(r — R) — (t — u){r' - R)} 2 
Ul ~l (t-u)(s-u) 2 
Vi = « . [(s-u)(r-R)-(t-u)(r'-R)] 2 
s (t — u) 2 (s — u) 

and an undetermined scaling function \l/ 3 . The G-functions express causality. 
Specifically, for the autocorrelator, i.e. r = r' this yields, with y = t/s 
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/g gfa-iy 1 (65) 

xa U(»-«)(l-«) 2 ' (»-«) a (l-«)J ( ' 

where if is an undetermined scaling function. A similar, but quite lengthy, expression 
can be derived for C 2 and depends on x 2 and ix [9J. Since the critical bcp is described 
by a free field-theory, one can expect from table [3] that x = x = d/2 and hence for the 
composite fields 

3 

x 2 = d, ix = ~d (66) 

Consequently, the autocorrelator takes the general form 

C(t, s) = s 1 - d / 2 g 1 (t/s) + s 2 ~ d g 2 (t/s) (67) 

For d larger than the lower critical dimension d* = 2, the second term merely furnishes 
a finite-time correction. On the other hand, for d < c?* = 2, it would be the dominant 
one and we can only achieve agreement by discarding the scaling function g 2 . It remains 
to be seen that gi is compatible with the exact result given in table HI 

This can be achieved by choosing in eq. fl63l) [104 J 

* 3 (u 1 ,v 1 )=e(—--) (68) 



Ui Vi 

where H remains an arbitrary function. Then 

d(t,s)= -ixpos-* +1 - x -^(t/s-l)h 



x J\e [(t/s -B)(i- o)]i-fr<k ( t/s t J s 1 _ 1 26 ) (69) 

where the function (fix is defined by 

0iH= fdRexp(-—R 2 )z(R 2 ) (70) 



Now the result for the bcp in table H] is recovered if one chooses [104[ [9] <pi(w) = 
(j) 0tC w~ 1 ~ a . This form for <fi\{w) guarantees that the three-point response function 
(4>(r, t)(j)(r, s)<p 2 {r\ u)) is nonsingular for t = s0 

jj We remark that for 2 < d < 4, the same form of the autocorrelation function is also found in the 
critical voter- model 371. 
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3.2. Bosonic pair- contact process 

The construction of the action follows standard lines (78]. Making the same shift eq. <^4 
as before, the action becomes J[0, 0] = Jq[4>, 4>] + Jb[4>, <f>] where the 'deterministic' part 
now reads 

J o [0, 0] : = Jdr J dt 4>(2Md t - V 2 )0 - «0 2 2 • (71) 



The remaining part is the noise-term which reads 
J b [<f>,$ = I dR I du 



2x2 - 2ap (f) 2 (j) - jj(f) 3 (f) 2 - 2/ip o 3 - " 2a3 



(72) 

The discussion of the Schrodinger- or, more precisely, the ageing-invariance of Jo can 
no longer use the representations we considered so far, since the equation of motion 
associated to Jo is non-linear, viz. 

2Md t <f){x, t) = V 2 <p{x, t) - g<p 2 {x, t)$(x, t) (73) 

While for a constant g the well-known symmetries of this equation are those encountered 
before, it was pointed out recently that g rather should be considered as a dimensionful 
quantity and hence should transform under local scale-transformations as well [119] . 
The requires an extension of the generators of agc d which do contain a dimensionful 
coupling g [9|. Then it can be shown that the Bargman supers election rules [2$) still 
apply and the response function of the noiseless theory now reads [9] 

-1(X!-X 2 ) 



Mt,s;r,r') = (t-s)-^ + ^ Q 



M ( r -r') 2 \ ~ ft t-s g , 



2 t-s J \s g 1 /y , {t-s)y^ 

with an undetermined scaling function \l/ 2 . In these calculation, we have assumed for 
technical simplicity that each field tpt has a coupling and only at the end, we let 
gi = . . . = g n . = g. This form is clearly consistent with the results for the bpcp listed 
in table HI for both a < a c and a = a c , if we identify 

x := X\ = x 2 = a + 1 = ^ , ^2 = const. (75) 

In distinction with the bosonic contact process, the symmetries of the noiseless part So 
do not fix the response function completely but leave a certain degree of flexibility in 
form of the scaling function ^2- 

As before, averages can be reduced to averages within the 'deterministic' theory 
only. Technically, calculations become a little longer for the bpcp, since because of the 
structure of J, several composite field must be defined. We refer to [9] for the details and 
merely quote here the results. First, the response function does not depend explicitly 
on the noise, viz. 



R(t,s;r,r f ) = R (t,s;r,r') (76) 
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Second, the results of table @] for fc(y) can be reproduced from the single term 

Ci(t, s) = ap\ JdR Jdu ($(t, r)<j>(s, r)0 2 (w, R)^ (77) 
The required three-point function now reads 



0(t, r)0(s, r')(j) 2 (u, R) ) = (t - s) x ~^ X2 (t - u)-**' (s - u)" 



-■>-'2 



x exp 



M(r- R) 2 M 



R) 



t — u 



u 



with 



Ml 



V i = - 



01 
S2 



u [(s - u )(r-R)-(t- u)(r' - R)} 2 
1 (t -u)(s- u) 2 

u [(s — u)(r — R) — (t — u){r' - R)] 2 
s 
1 



(t — u) 2 (s — u) 
a 1 ^ „ 1 a 1 ^ 



s 2 (t-u) 2 ' 
1 



02 



S 2 



U 



1 

+ - 



t — u u 

Next, we choose the following realization for 



^3(ui,v h 0i,0 2 ,03) 



^ 3 {u 1 ,v 1 ,0 1 ,0 2 ,0 3 ) (78) 

(79) 

(80) 

(81) 
(82) 



03 = a l ' y s 2 



i 1 - 









Wi~ V02)V0; 

03-VM~3 



(a-b) 



(83) 



where the scaling function S was already encountered in eq. (1681) for the bosonic contact 
process. We now have to distinguish the two different cases a < «c and a = ac- For 
the first case a < ac, we have a — b = so that the last factor in flHBl disappears and 
we simply return to the expressions already found for the bosonic contact process, in 
agreement with the known exact results. However, at the multicritical point a = ac 
we have a — b 7^ and the last factor becomes important. We point out that only the 
presence or absence of this factor distinguishes the cases a < ac and a = ac- 

Inserting the values of the /3i,2,3 we finally obtain for the autocorrelation function 

d{t, s) = s-\y - f d9 [(t/s - 0)(1 - 6)} a - b 

Jo 

ft/s + i -26\ r mi si - n ~\ a ~ h 

X 01 ' 



t/s 



6(t/s - 1) 



(t/s-6)(l-6) 



where we have identified 

x 2 = 2(& — a) + d 

C\(t,s) reduces to the expression ( |53l if we choose the same <pi(w) = 4>o, c w 
before. Hence all scaling functions for the bpcp are reproduced correctly. 



-l-a 



4) 



55) 



as 
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4. Growth models 

As a further illustration on a different class of system we now describe recent work by 
Rothlein, Baumann and Pleimling [110] on kinetic growth models. As a simple model 
for ballistic deposition consider the Family model [H]. A particle is randomly dropped 
onto the sites of a lattice. However, before it is fixed, the particle explores the sites 
around the one it arrived at (typically the nearest neighbours) and fixes itself at the 
lattice site with the lowest height. One obtains in this way a growing surface which may 
be described in terms of a height variable h(t, r). Clearly, since there is only irreversible 
deposition, the system will never arrive at an equilibrium state. 

Some of the simplest continuum descriptions of these phenomena can be cast into 
stochastic linear equations for the height variable. For the sake of simplicity, we shall 
always work in the frame co-moving with the mean surface height. For example, if the 
deposition is purely diffusive and without mass conservation, the simplest model is the 
well-known Edwards- Wilkinson (ew) model [ID] 

d t h(t, r) = DX7 2 h(t, r) + rj(t, r) (86) 

but on the other side, if mass conservation must be taken into account, one might rather 
consider the Mullins- Herring (mh) model, see [127] 

dth(t, r) = -Diy 2 fh{t, r) + r)(t, r) (87) 

Following [110J, the following types of gaussian noise with vanishing first moment 
(rj(t, r)) = will be considered: 

(a) non-conserved, short-ranged (f](t, r)r)(s, r')) = 2D5(r — r')8(t — s). 

(b) non-conserved, long-range (r)(t,r)r)(s,r')) = 2D\r — r'\ 2p ~ d 5(t — s) and < p < d/2. 

(c) conserved, short-ranged (r)(t, r)r)(s, r')) = —2DV 2 5{r — r')S(t — s). 

Then the following models were studied in [110] : 

(i) EWl: eq. fl86|) with the non-conserved noise (a). 

(ii) EW2: eq. ( 1561 with the non-conserved, long-ranged noise (b). 

(iii) MHl: eq. (1571) with the non-conserved noise (a). 

(iv) MH2: eq. (1571) with the non-conserved, long-ranged noise (b). 

(v) MHc: eq. (1571) with the conserved noise (c) [14] . 

In the models EWl and MHc the noise is in agreement with detailed balance while for 
the other models it is not. Solving the linear equations (156|) and (157}) is straightforward. 
The two-time correlation and response functions are seen to obey the same kind of 
scaling behaviour as for the non-equilibrium models considered before. In table [3j the 
values of the exponents are listed and the scaling functions for the autoresponse and 
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the autocorrelation are included in table H] (the models bcp and EWl lead to identical 
results and are not listed separately). The result quoted for the MHc model is only valid 
for d > 2 as stated; for d = 2 the scaling function becomes fc(y) = 2D ln[(y — l)/(y+ 1)] 
[TJ] • Detailed simulations show that the correlation and response functions of the Family 
model [B] and of a variant of it are perfectly described by the EWl model and hence 
should be in the same universality class jllOj . 

We shall see shortly that all results about the scaling functions as listed in table H] 
can be understood from an extension of dynamical scaling to local scale-invariance, for 
both the EW models, where the dynamical exponent z = 2, as well as for the MH models, 
where z = 4. In this context, it is instructive to consider the space-time responses 
R(t, s; r) as well. In general, one finds the structure 



R{t,s;r) = R(t,s)$(\r\(t- s) 



-X/z\ 



where R(t, s) = R(t, s; 0) is the autores pon se function. If z — 2, one has a simple 
exponential form = exp(— A4u) [56_|jjj while for z — 4, the form of $(w) is more 
complicated. For example, in the MHc model with conserved noise [14] 



2'4'256y rfr(f + |)V 256 V 2 '4 2' 256 



(89) 



where $o is a known constant and 0-^2 is a hypergeometric function. Very similar 
expressions have been derived in the other MH models, see |110j for details. 

Understanding the form of the spatio-temporal response allows for an explicit test 
of one of the important ingredients of the theory of local scale-invariance, namely Galilei- 
invariance for z = 2 and its generalization if z 7^ 2. Since in the EW models one has 
z = 2, the calculation of R(t, s; r) is a direct extension of the discussion presented in 
section 1. Since this has been discussed in detail in the litterature [561 EH] we shall not 
repeat it here and rather concentrate on the case z = 4, following [581 1110] EH]- 

Consider the dynamical symmetries of the 'Schrodinger operator' 

S 4 :=-Xd t + ^(V 2 r ) 2 . (90) 

which will become related to the deterministic part of the Mullins-Herring equation. 
As before, we ask if dynamical scaling can be extended to a larger set of local scale- 
transformations, given standard dynamical scaling and spatial translation-invar iance. 
Specializing the construction of [SS] to z = 4, these generators read as follows, with the 



shorthands r ■ d r := J2k=i r k dr k ^l ■= J2k=l d r k and r<2 '■= J2k=l 



ki 



ff For detailed quantitative tests in Ising and Potts models, see [5T1 155] . 
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X-! := -d t 

X Q := -tdt - -r ■ d r - - 

X, := -t 2 d t -^t-Xr 2 (Vl)- 1 -^tr.d r 

+ 4 7 (r • d r ) (Vj) ~ 2 + 2 7 (rf - 4) (V 2 r ) ~ 2 (91) 

jib J) := nd r . — Tjd ri ; where 1 < i < j < d 

yW - _ f) 
I -i/4 ~ u n 

1$ = -^ ! -4Ar 4 (V^)- 1 + 8 7 a ri (V^)- 2 

where x is the scaling dimension of the fields on which these generators act and 7, A are 
further field-dependent parameters. Here, the generators X± 1)0 correspond to projective 
changes in the time t *—>■ (at + /3)/(t£ + 5) with a5 — (3^ = 1, the generators Y^-im are 
space-translations, generalized Galilei-transformations and so on and R^'^ are spatial 
rotations. Here we use the following properties d^d@ = d® +13 and [<9",r] = ad" -1 , which 
can be proven for fractional derivatives with extra distributional terms [581 Appendix A] 
which in turn are closely related to fractional derivatives as defined in |48j . Furthermore, 
(V^)~ 2 = (V 2 ) -1 • (V 2 ) -1 and the operator (V 2 )" 1 is defined, e.g., for d = 2, by formal 
expansion [1101 [1] 



n=0 

This implies the commutator [(V^) n ,rj] = n<9 rj (V 2 ) n_1 for all n e Z. 

That the generators eq. fl9Tl) indeed describe dynamical symmetries of the 
'Schrodinger operator' f[9"Ul) now follows from the commutators 



[5 4 ,ri? /4 ] = [s,,y^\ = [«s 4 ,^] = o , [s,,x \ = -<s 4 (93) 

This means that for a solution of the 'Schrodinger equation' S^cf) = the transformed 
function Xcf) is again solution of the 'Schrodinger equation'. Finally 

[<S 4 , X x ] = -2tS 4 + Ux - i - 1 + ?1\ (94) 

hence a dynamical symmetry is obtained if the field <p has the scaling dimension 
d 2 7 

x = 2 + 1 " T (95) 

Generalizing from conformal or Schrodinger-invariance, quasiprimary fields transform 
covariantly under the generators ( 19TI) and in particular the response function will satisfy 
the conditions X\R = X R = Y_y A R = R^R = (the other conditions then follow 
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from the Jacobi identities) and is now characterized by its scaling dimension Xj and the 
further parameters 7$, Aj. In calculating the response function R = (<fi<f)) this leads to 
the conditions A = —A and 7 = —7 whereas the scaling function from eq. ( 1881) can 
be found by solving the differential equation 

d u (Jr~A K" 1 ^)) + 4Am {Ji-A " 16 ^«j = (96) 

where 




Solving this via series expansion techniques [58l 11101 [14] and checking carefully that all 
independent solutions are taken into account, one can indeed recover the explicit result 
( |89|) for the MHc model as a special case [14] and, similarly, also for the MHl and mh2 
models [1 10J . 

Lastly, since the MH equation is linear, there is a natural Wick theorem which allows 
to go over from the stochastic Langevin equation to the deterministic equation in quite 
an analogous way as previously for z = 2 [1101 [14] . The extension of the technique to 
non-linear cases and/or to z 7^ 2,4 is work in progress and will be reported elsewhere 
[15] . In a similar way one may also check that the correlation functions agree with LSI. 

The MH models considered in this section are, together with the critical spherical 
model with a conserved order-parameter [TJ] , the first analytically solved examples with 
z 7^ 2 where local scale-invariance could be fully confirmed. These examples make it in 
particular clear that the height of the surface in growth processes is a natural candidate 
for being described by a quasiprimary scaling operator of local scale-invariance. 



5. Conclusions 

In this survey, we have reviewed to what extent one may expect that a phenomeno logical 
description, which has been successfully applied to describe the ageing of magnetic 
systems relaxing towards equilibrium steady-states, may be extended to more general 
models where the stationary states are no longer part of an equilibrium statistical 
ensemble. This situation frequently arises in chemical kinetics, see [H [6j [311 I124[ 
125], and the consideration of such systems is of interest in studies of ageing in 
chemical/biological systems where already the intrinsic microscopic dynamics and/or 
constraints does not admit relaxation towards thermal equilibrium. It is clear that the 
study of ageing phenomena without detailed balance still stands at its very beginning 
and many open questions remain. In particular, the few models reviewed here certainly 
do not exhaust all possibilities for ageing behaviour without detailed balance but should 
be rather seen as case studies whose results might suggest further research problems. 
This review has already served a useful purpose if it encourages people to explore more 
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systematically the properties of two-time observables of non-equilibrium systems, e.g. 
react ion- diffusion problems or biologically motivated models. 

Specifically, the following points should be noted: 

(i) the generic scaling form (|28|) was seen to be satisfied in all models considered. 
However, the exponent relation a = b, known to hold for critical systems with 
detailed balance, is no longer valid in general, see table [31 For different universality 
classes, the relation a = b is either maintained or else broken in different ways. 
This means that there is no obvious and general analogue of an universal limit 
fluctuation-dissipation ratio (for magnets one sometimes tries to relate this 
to a non-equilibrium temperature, see e.g. [33]) even if such an analogy may be 
defined for certain subclasses. 

(ii) for the uncorrelated initial states which have been considered so far, one observes 
that the autocorrelation and autoresponse exponents agree Ac* = A#. It would 
be interesting to see if spatial or temporal disorder in the rates may change that 
conclusion, as it apparently happens in diluted magnets [114[ [72]. 

(iii) one of our main question with respect to ageing systems has been if dynamical 
scaling permits an extension to a larger group of local scale-transformation |58j . 
It has turned out that from the point of view of LSI the responses are the most 
easy quantities to study. The long list of examples, see tables [T] and [2], where the 
two-time autoresponse function R(t, s) was concluded to be in agreement with LSI 
is clear evidence that LSI is indeed a successful phenomenological scheme, and this 
for values of the dynamical exponent z which are often far from z = 2 characterizing 
simple diffusive motion. But this ansatz remains to be proven, especially for z ^ 2, 
e.g. from some underlying stochastic Langevin equation. The exactly solvable 
examples we have treated suggest that the idea of splitting the Langevin equation in 
a 'deterministic' part with possible non-trivial dynamical symmetries and a 'noise' 
part which breaks those can be taken over from magnets to more general reaction- 
diffusion type system, although the noise terms can be considerably more general. 
The examples studied here also suggest that the basic physical variables of these 
models, such as the particle-density or the height of the surface, should be directly 
relatable to the quasiprimary scaling operators of LSI. 

It appears to us that it should be promising to investigate more systematically the 
foundations and consequences of a hitherto unsuspected non-trivial dynamical symmetry 
in scale-invariant non-equilibrium dynamics. 
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